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ABSTRACT 

Rate coefficients for rotational transitions in H 2 induced by H 2 impact are 
presented. Extensive quantum mechanical coupled-channel calculations based 
on a recently published ^2)2 potential energy surface were performed. The 
potential energy surface used here is presumed to be more reliable than surfaces 
used in previous work. Rotational transition cross sections with initial levels J < 
8 were computed for collision energies ranging between 10~ 4 and 2.5 eV, and the 
corresponding rate coefficients were calculated for the temperature range 2 < T 
< 10,000 K. In general, agreement with earlier calculations, which were limited 
to 100-6000 K, is good though discrepancies are found at the lowest and highest 
temperatures. Low-density-limit cooling functions due to para- and ortho-H2 
collisions are obtained from the collisional rate coefficients. Implications of the 
new results for non-thermal H 2 rotational distributions in molecular regions are 
also investigated. 

Subject headings: molecular processes — - molecular data — ISM: molecules 

1. INTRODUCTION 

Collision-induced energy transfer involving H 2 molecules plays an important role in many 
areas of astrophysics since hydrogen is the most abundant molecule in the interstellar medium 
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(ISM), dominating the mass of both dense and translucent molecular clouds. Because H2 is 
a homo nuclear molecule and lacks a dipole moment, its infrared emission is weak and it is 
difficult to observe directly. Furthermore, because of the low mass of the hydrogen atom, 
the rotational levels of molecular hydrogen are widely spaced and require relatively high 
temperatures, T > 100 K, to excite even the lowest levels appreciably via thermal collisions. 
Nevertheless, with the current sensitivity of detectors deployed on space-based observatories, 
such as the Spitzer Space Telescope and the former Infrared Space Observatory (ISO), the 
detection and mapping of the pur e rotational lines of H 2 in a significant sample of galaxies is 



now possible (INeufeld et al 



from the ultra yiolet (IMever et al 



near infrared (IBurton et al. 



2006 



1998 



2001 



Fuente et al.lll998l). These observations, as well a s others 



McCartney et al 



Boisse et al. 2005; France et a 



1999 



Allers et al 



20051 . 120071) to the 



20051 ) probe shock- 



induced heating and cooling and UV and x-ray irradiation of gas in the ISM of galactic and 
extragalactic sources and star-forming regions. 

Being the simplest diatom-diatom collision system and one in which all constituents 
are quantum mechanically indistinguishable, early work, notably by Green and coworkers, 
carried out a succession of theoretical studi e s on rotational re laxation in H 2 +H 2 collisions 
(jGreenlll975l ; iRamaswamy et al.lll977l . 1 19781 ; iGreen et al.lll978l ). The calculation of the re- 
laxation cross sections involve the quantum mechanical scattering of the heavy particles on 
a potential energy surface (PES) which is computed using quantum chemistry techniques, 
accounting for the electronic motion for all relevant nuclear configurations of the constituent 
colliders. The availability of more realistic (H 2 )2 potential surfaces has enabled consider- 
able progress in obtaining a more reliable set of state-to-state collisional cro ss sections and 
rate coefficients. To d ate, t he most comprehensive study has been made by iFlowerl (119981 ) 
and iFlower fc Rouefil (119981 ) . Within the rigi d-rotor approxim ation, Flower & Roueff have 
made use of the (H 2 ) 2 interaction potential of ISchwenkd (119881 ) in a quantal coupled-channel 
method to determine the rate coefficients for rotational transitions in H 2 + H 2 collisions. 
Rotational levels J < 8 and kinetic temperatures T < 1000 K were employed and the results 
were compared with those reported by Danby, Flower, & Monteiro (1987) based on an older 
potential surface. 

Accurate determinations of rate coefficients for state-to-state transitions in small molec- 
ular systems such as the H 2 molecule require a quantum mechanical description of the scat- 
tering process. A full quantum calculation of rovibrational transitions in the H 2 -H 2 system 
is computationally challenging. So far, only a few studies have been reported that include 
both rotational and vibrational degrees of freedom of the H 2 molecules. Pogrebnya and Clary 
(2002) reported vibrational relaxation of H 2 in collisions with H 2 using a coupled-states ap- 
proximation implemented within a time-independent quantum mechanical approach. Panda 
et al. (2007) and Otto et al. (2008) employed a time-dependent quantum mechanical ap- 
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proach and reported rotational transition cross sections in ortho-H2+para-H2 and para-H2 
collisions. More recently, Quemener et al. (2008) reported a full-dimensional quantum scat- 
tering calculation of rotational and vibrational transitions in the H 2 -H 2 system that does not 
involve any dynamics approximation. All these calculations adopted the six-dimensional H 4 
PES developed by Boothroyd et al. (2002). While this PES is useful for benchmarking full- 
dimensional quant um scattering codes for the H2-H2 system, co mparisons with experimental 
data have shown (ILee et al.l 120061 ; lOtto. Gatti. fc Meyerl 120081 ) that the PES predicts rota- 
tional transition rate coefficients between the J = and J = 2 levels that are too small and 
that it is not appropriate for accurate determination of rotational transition rate constants 
in the H 2 -H 2 system. 

Since full-dimensional quantum calculations of four-atom systems are computationally 
demanding, especially for excited rotational levels, the rigid-rotor approximation is often 
employed. The adequacy of the rigid-rotor approximation for pure rotational transitions in 
H 2 -H 2 collisions was recently demonstrated by Otto et al. (2008) who reported cross sec- 
tions for rotational excitations in para-H 2 using the rigid-rotor model and a full- dimensional 
quantum calculation. While these calculations have been performed using the BMKP PES, 
we believe that the conclusions are valid for the H 2 -H 2 system in general. 

The accuracy of collisional data for astrophysical applications is limited by the uncer- 
tainty in the PES and the approximations employed in the dynamics calculations. It is a 
non-trivial task to judge the impact of the uncertainties in the collisional data, associated 
with the choice of the PES. Accurate determination of the potential energy in the interaction 
region to the requisite level of accuracy of ~ 10~ 3 hartrees (or ~10 K) remains a challenge, as 
mentioned above. This is especially important at low temperatures where small uncertainties 
in the interaction potential translate into larger errors in the cross sections and rate con- 
stants. A further difficulty arises from fitting a limited number of ab initio potential energy 
points to obtain a surface over all configurations which are sampled in the scattering calcula- 
tions. Unphysical behavior in the fitted PES can result, particularly when bridging explicitly 
calculated energies at intermediate internuclear separations to long-range asymptotic forms. 



Recently, a new and improved (H 2 ) 2 rigid- rotor PES was computed by lDiep fc Johnson 



( 120001 ). The accuracy of this PES for rotational transitions in H 2 has recently been demon- 
strated by comparing computed rate coefficients for J = — > 2 rotational excitation against 
experimental results (IMate et al.ll2005l ; ILee et al.ll2006l ; lOtto. Gatti. &: Meyerl I2008I ). Thus, 
we believe that the new PES by Diep & Johnson (2000) could be employed to provide re- 
liable values of rate constants for rotational transitions in the H 2 +H 2 system. Since the 
H 2 +H 2 collision system is of astrophysical significance, it is crucial to establish whether the 
cross sections and rate coefficients for rotational energy transfer based on this PES are in 
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agreement with the earlier data. Therefore, we adopted this PES in the present quantum 
mechanical close-coupling calculations to obtain the rotational transition cross sections with 
initial levels J < 8 for collision energies ranging from 10~ 4 to 2.5 eV. The corresponding 
rate coefficients were computed for the temperat ure range of 2 < T < 10,000 K. The presen t 
rate coefficients are compared with the results of iFlowerl (119981 ) and lFlower fc Roueffl (119981 ). 
which are to the best of our knowledge, the only comprehensive calculations available for 
rotational de-excitation rate coefficients. We also present H 2 -H 2 cooling functions in the 
low-density limit and test the new rate coefficients in simulations of UV irradiated molecular 
gas. Atomic units are used throughout unless otherwise specified, i.e., distances are in bohrs 
(a ) and energy in hartrees (E h ). Recall that 1 a = 0.529177 A, while 1 E h = 27.2114 eV 
= 219474.635156 cm" 1 = 627.51 kcal/mole. 



2. COMPUTATIONAL DETAILS 

We carried out quantal coupled-channel calculations for collisions of H 2 with H 2 us- 
ing a fixed bond length of 1.449 a Q (0.7668 A). The rigid-rotor H 2 target and projectile 
energy levels were calculated using a rotational constant of B = 60.853 cm -1 . To solve 
the coupled-channel equations, we used the hybrid modified log-derivative Airy propaga- 
tor (Alexander & Manolo poulos 1987) in the gen eral purpose non-reactive scattering code 



MOLSCAT developed by iHutson fc Greenl ()1994j). The log- derivative matrix is propagated 



to large intermolecular separations where the numerical results are matched to the known 
asymptotic solutions to extract the physical scattering matrix. This procedure is carried out 
for each partial wave. We have checked that the total integral cross sections are converged 
with respect to the number of partial waves, as well as the asymptotic matching radius (e.g., 
R = 40 a ) for all channels included in the calculations. 

In addition to the partial wave convergence, we have checked that the results are op- 
timized with respect to other parameters that enter into the scattering calculations. In 
particular, the parameters used for the analytical expression for the Diep & Johnson PES, 

v(r, 0!, e 2 , 12 ) = v hM {R)G hM {e u e 2 , 0i 2 ), (i) 

h,h,l 

where Vi lt i 2t i(R) are radial expansion coefficients and Gi lt i 2 j{0\, 6* 2 , 0i 2 ) are bispherical har- 
monics. The angles 9i, 6* 2 denote the two in-plane angles and i2 is the relative torsional 
angle. We used 10 quadrature points each for integration along angular coordinates 9i, 02, 
and <pi2- From the Diep & Johnson fit to their numerical PES, the expansion coefficients 
Vi u i 2t i(R) are determined. They noted that only the Vo,o,o(-R)> ^0,2 (-R), Vo, 2)2 (-R) and V^^i?) 
terms make significant contribution to the potential. 
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In Table [p we show a comparison between the PES adopted by Flower & Roueff, the 
Schwenkd ( 119881 ) potential, and the Diep & Johnson (H 2 )2 PES for linear [9i, 9 2 , </>i2] = [0, 
0, 0] and parallel [9\, 9 2 , (P12]— [it /2,tt/2,0] configurations, respectively. Although the two 



potentials show close agreement for both configurations, the analytic PES of lSchwenke Jl988h 
for (112)2 is known to be accurate only for pairs of hydrogen molecules with intermolecular 
separations not greater than ~ 5 a Q . The small discrepancy may be attributed to the 
improved accuracy in the Diep & Johnson PES which incorporates the correct long-range 
behavior. 

The present calculations consider the hydrogen molecules to be indistinguishable. Sym- 
metric basis sets were used for para-para and ortho-ortho collisions, whereas, for para-ortho 
collisions, asymmetric basis sets were chosen based on the order of the energy spectrum. 
The scattering cross sections for rotational transitions were computed for collision energies 
ranging between 10~ 4 and 2.5 eV. Four de-excitation cases with AJ\ = —2 were considered 
(i.e., para-para and para-ortho; ortho-para and ortho-ortho). Note that a sufficiently large, 
but truncated basis set has been used to optimize the computation time and minimize the 
loss of numerical accuracy (i.e., to < 5%) in th e scattering cross section calculations. Further 
details about the calculations can be found in lLee et al.l (120061 ). 



Rate coefficients for state-to-state rotational transitions were obtained by averaging the 
appropriate cross sections over a Boltzmann distribution of the relative kinetic energy E k of 
the H 2 molecules at a given temperature T 



k 



(T) 



G 



dE h a Jx 



J2- 



*j[j' 2 (E k )E k e 



(2) 



where G 



ft = {ksT) 1 with ks being the Boltzmann constant, and we have 



adopt ed the definition of the cross section for indistinguishable particles as given by iGreen 
(119751). This defin i tion d iffers from that a dopted by Flower fc Roueff which was given 
by Zarur fc Rabitz ( 1974 ). As discussed by Danby. Flower, fc Monteirol ( 1987 ). for single 
rotational excitation transitions, the cross section defined by I Green £973) must be divided 
by two if J\ = J2 or J[ = J' 2 to prevent double counting in the determination of the 
production rate of J[ or J' 2 - The factor (1 + 5j 1 j 2 )(1 + <5j(j^) in the denominator of eq. (j2j) 
accounts for it. 
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RESULTS 



3.1. Comparison of Rate Coefficients 

Fig. 1 shows a comparison betw een the present theoretical results and the experimental 
rate coefficients of lMate et al.l (120051 ) for the (0,0)— »(2,0) excitation transition between 50 and 
300 K, where we have used the notation (Ji, J2). The experimental and current theoretical 
rate coefficients are in good ag reement as wa s shown previously in Lee et al. (2006). The 
theoretical results obtained by iFlowerl (119981 ). which used the older PES of Schwenke, also 
shows good agreement with the measured and new calculated rate coefficients. 

Figures 2(a) and 2(b) show the rotational de-excitation rate coefficients AJi = —2, 
AJ2 = for para-para and ortho-ortho collisions, respectively. For the para-para case, the 
current rate c oeffic ients for most transition s are found to be similar to those obtained by 
Flower! (119981 ) and iFlower fc Rouen (119981 ). but for the important (2,0)— >(0,0) transition, 
Flower & Roueff's result is ~10% smaller for tem peratures less tha n 1000 KLj. On the other 
hand, (8,0)— ^(6,0) transition rate coefficients from lFlowerl (119981 ) and lFlower &: Roueffl (119981 ) 



are larger compared to the current results and the agreement appears to deteriorate as the 
temperature increases above ~400 K. The (4,0)— >(2,0) transition, however, is in excellent 
agreement, but only below ~1000 K. Likewise, for the ortho-ortho transitions, the rate 
coefficients of Flower & Roueff are similar to the current results, except at the highest 
temperatures. 

Figures 3(a) and 3(b) show the rotational de-excitation rate coefficients for para-ortho 
collisions. Reasonable agreement was observed, except for the (8,1)— »(6,1), where the results 
of Flower & Roueff tend to be significantly larger than the present calculations. This be- 
havior is essentially identical to that seen for the (8,0)— >(6,0) transition shown in Fig. 2(a). 
Discrepancies are also noted for all of the transitions at the highest temperatures. 



3.2. Rate Coefficient Data and Fitting 

Numerical values for the rate coefficients are given in Table 2. We have also fit- 
ted the computed rate coefficients to an analytical form similar to that introduced by 



We note that the rate coefficients from 



Flower (1998) are limited to temperatures between 100 and 



1000 K, while those given in blower fc Rouefj (|l998h extend from 1000 to 4500 K, but only for para-H 2 
colliders. We supplement these in Figs. 2 and 3, and in the discussion throughout, with additional data 
given between 1500 and 6000 K on David Flower's website at ccp7.dur.ac.uk/cooling_by_h2/nodel.html. 
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Le Bourlot et al.l (119991 ). The adopted form used here is 

kj lJ2 ^ J[4 (T) = W (*+b/ti+c/ti) + 10 [e+(//t 2 )"] (3) 

where t = T/1000 + 5t, t x = d x T/1000 + 8t, and t 2 = g x T/1000 + 5i, with 6t = 1. The 
fit coefficients a through /i, a Fortran routine to return the rate coefficient k following input 
of the temperature T and the set of initial and final quantum numbers, the numerical data 
for the rate coefficients, and the numerical data for the cross sections can all be obtained at 
the website: www.physast.uga.edu/ugamop/. The fits are valid between 2 and 10,000 K and 
were performed for exothermic, AJ X = —2, transitions. Rate coefficients for endothermic 
transitions can be obtained by detailed balance which is implemented in the Fortran routine. 

Because the accuracy of the PESs for values of li, l 2 , and I > 4 (see eq. 1) is uncertain, 
we did not explicitly compute the rate coefficients for AJ\ < —2 transitions, except for two 
AJi = —4 cases: (4,0)— >(0,0) and (5,1) — ^(1,1). The rate coefficients for these processes are 
typically a factor of ~50 smaller than the AJi = —2 rate coefficients, for the same initial 
state. Transitions with A J\ < —4 are expected to have even smaller rate coefficients. While 
these two rate coefficients have been fitted, we estimate all other A J% < —2 transitions with 



a version of the Modified Exponential Energy Gap Law (MEGL) (ISteinfeld et al.l Il99ll ). 
With the MEGL, the rate coefficients for AJi < —2 transitions are obtained by scaling the 
AJi = —2 rate coefficients with the following relation 

t (T\ u ^ 2(^-2) + ! Qfrr ^ /rr ^ l+aA/T f 

k Jx j^j [Jk {T) = fc JlJWl _ 2i4 (r) — — — — exp(-/3(T)A/T) ^ + ( 4 ) 

where A = [E(J\ — 2) — E(J[)]/kB is the difference in rotational binding energies in K and 
a = 2 is usually adopted. We find that the (4,0)— >(0,0) transition is fit with 7 = 2/3 and 
f3{T)/T = 1/300, the latter being nearly temperature independent. 



4. Astrophysical Applications 

In order to investigate the effect of the new rate coefficients in astrophysical environ- 
ments, we show in Figure 4 low-density-limit cooling functions given by 

A(T) 

where A(T) is the cooling rate in erg cm" 3 s -1 , % 2 (J) is the number density of H 2 in 
rotational state J, and huj^ is the emitted photon energy. Here we consider transitions 
among the lowest four rotational states and plot the cooling rate due to para-H2 collisions 
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on H 2 (Ji = 0, 1) and the total cooling rate for a statistical dis tribution of ort h o- an d para- 
H 2 , i.e. ortho:para=3:l, following the prescription given by iGlover fc Abell (|20081) . The 



new cooling functions are compared to those obtained by IGlover fc Abell (120081 ). based on 
the H 2 -H 2 rate coefficients of Flower & Roueff, which are indicated to be valid between 
100 < T < 6000 K. The two sets of cooling functions are generally in agreement to within 
~20% for T < 2000 K including cooling due to ortho-H2 collisions (not shown). The level 
of agreement reflects differences in the collisional rate coefficients show n in Figs. 2 and 3. 
However, above ~2000 K the cooling functions of IGlover fc Abell (120081 ) are seen to depart 
significantly from the current results which is due to the divergence in the two sets of rate 
coefficients at the higher temperatures. 

In a further test, we include the full set of H 2 -H 2 rate coefficients in a model of a 



mole cular photodissociation region using the spectral si mulation code Clou dy (IFerland et al. 



19981 ). Model V3 of the PDR-Code Comparison Study (iRollig et all 120071 ) is adopted where 
the total density is taken to be 10 5,5 cm -3 and t he inten s ity of the incident far ultraviolet 
(FUV) radiation field is 10 times the standard iDraind (119781) FUV field. Other model 
parameters of the molecular region are given in IRollig et al.l (120071 ) . A nearly complete 
description of the electronic, vibrational, and rot ational radi a tive a nd collisional properties 
of the H 2 molecule was included as described by IShaw et al.l (120051 ) . Figure 5a gives the H 



and H 2 number densities and the temperature structure of the cloud as a function of depth. 

In Figure 5b, the effect of the new H 2 -H 2 rate coefficients are shown by plotting the 
ratio of the rotational populations n(v = 0, J) determined from two sep arate models. I n 



the n umerator, the standard set of H 2 -H 2 rate coefficient fits given in iLe Bourlot et al. 



( 119991 ) are adopted, while in the denominator the current H 2 -H 2 rate coefficients are used. 
H 2 -H and H?-He collisions ar e also incl uded and t aken from the recent calculations of 
Wrathmall. Gusdorf. fc Flowerl (120071 ) and iLee et al.l (120081 ) and used in both models. Since 
the rate coefficients from the two calculations are in general agreement, most levels do not 
experience a significant change in population. However, differences are seen for cloud depths 
greater than ~ 10 10 cm which corresponds to an increase in the total H 2 abundance and 
decreases in the H abundance and gas temperature. The abundance changes maximize the 
effect of H 2 colliders reducing the importance of H collisions, though He collisions are still im- 
portant. As such, while the discrepancies in the H 2 -H 2 rate coefficients are maximum above 
~1000 K, their effects at these high temperatures will be masked by the more important roles 
of H and He collisions. However, as the temperature decreases below 100 K and H 2 collisions 



becom e more important, the discrepancies are related to the fact that the ILe Bourlot et al. 



( 119991 ) fits were not intended to be valid below 100 K. The largest differences are for the 
J = 3 and J = 4 populations, ~50% and ~20%, respectively. These differences are related 
to the increase in the rate coefficients to populate the levels as shown in Figs. 3a and 2a, re- 
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spectively, for temperatures less than 200 K. As the temperature falls further, below ~20 K, 
discrepancies in the J = 2 and J = 5 populations also become significant due to the fact 
that the current r ate c oefficients reach a minimum near this temperature, while the fits of 



Le Bourlot et al.l (119991 ) go to constant values. 



5. Summary 

We have performed an extensive quantum mechanical coupled-channel calculation, based 
on a recent and presumably more reliable (112)2 potential energy surface, to obtain rotational 
transition cross sections with initial levels J < 8. Collision energies ranging between 10~ 4 
and 2.5 eV were considered. We have computed the corresponding rate coefficients for the 
temper ature range fr om 2 to 10,000 K and have compared our results with the previous 



work of iFlowerl (119981 ) and lFlower fc Rouefil (119981 ). which were based on an older potential 



energy surface. We conclude that in some cases, the improvement made in the new potential 
energy surface and larger number of basis functions used in the present scattering calculations 
led to significant changes in the resulting state-to-state cross sections and rate coefficients. 
Implementation of the new H2-H2 rate coefficients in calculations of the low-density-limit 
cooling function and in the rotational level populations in a molecular photodissociation 
region result in some differences compared to the use of the rate coefficients from Flower & 
Roueff. Complete tabulations of rate coefficients are available electronically upon request, 
on the web at www.physast.uga.edu/ugamop, and from the ApJ electronic tables. 
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(NB), NASA grant NNG06GC94G (NB), NSF grant PHY-0554794 (RCF), and NSF grant 
AST-0607733 (PCS). We thank the NSF funded Institute for Theoretical Atomic, Molecular, 
and Optical Physics at the Harvard-Smithsonian Center for Astrophysics for travel support. 
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Temperature (K) 



Fig. 1. — Rate coefficients for (0,0)— >(2,0) rotational excitation in indistinguishable H 2 +H 2 
collisions as a fu n ction of temperature. The solid curve denotes the present resu lt using the 



Diep &: Johnson! (120001 ) PES. The filled squares are t he exper i ment al data from lMate et al. 



(120051 ) . The open circles are the theoretical results of iFlowerl (119981 ). 
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Fig. 2. — Temperature dependence of rotational state-resolved rate coefficients for H2+H2 
collis ions, (a) Par a-par a and (b) ortho-ortho c ollisions. Circles denote theoretical results 
from lFlower! (Il998h and iFlower fc Roueflj Jl998h . 
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Fig. 3. — Same as in Fig. 2, but for para-ortho collisions. 
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Temperature (K) 



Fig. 4. — Low-density-limit cooling func tion of due to H?+ rF. collisions. Thick lines show 
the present results. Thin lines are from I Glover fc Abel! ( 120081 ) based on the rate coefficients 
of Flower & Roueff. Solid, para-para; dash, para-ortho; and dot-dash, statistical ortho-para 
population. 
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Fig. 5. — (a) Total H and H 2 abundances and the gas temperature as a function of depth into 
a model molecular photodissociation region. See text for model parameters, (b) Effect of new 
H2-H2 rate coefficients on the H2 rotational populations in the model mole cular photodisso 



ciation region. Ratio of n-a 2 (v = 0, J = — 8) obtained with the standard iLe Bourlot et al. 



(119991 ) rate coefficient fits (based on Flower 1998) to those obtained by replacing the H 2 -H 2 
rate coefficients with the current calculations. See text for details. 
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Table 1: H 2 -H 2 interaction energy (in hartrees) for [6*i, 9 2 , 0i2] = [O, 0, 0] (upper entries) and 
02, <fii 2 ] = [ir/2, 7r/2, 0] (lower entries) configurations, respectively. 



^(Bohr) 


Flower (1998) 


Diep & Johnson (2000) 


3.0 


0.06546 


0.061703 


3.5 


0.02836 


0.027103 


4.0 


0.01168 


0.011378 


4.5 


0.004593 


0.004525 


5.0 


0.001687 


0.001660 


3.0 


0.03906 


0.038403 


3.5 


0.01679 


0.016568 


4.0 


0.006942 


0.006585 


4.5 


0.002761 


0.002343 


5.0 


0.001021 


0.000668 



Table 2: Rotational deexcitation rate coefficients as a function of temperature for H 2 (i>i = 
0, Ji) due to collisions of H 2 (i> 2 = 0, J 2 = or 1). 



Vi 


Ji 


v[ J[ v 2 J 2 v' 2 J' 2 T (K) 


k (cm 3 s *) 





2 


2.0 


8.365(-13) a 





2 


3.0 


7.416(-13) 





2 


4.0 


6.531(-13) 


The remainder of this table is 


available only on- 


-line as a machine readable table 





"The notation A(-B) refers to A x \Q- B . 



